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We study the dynamics of quantum excitations inside macromolecules which can undergo conformational 
transitions. In the first part of the paper, we use the path integral formalism to rigorously derive a set of coupled 
equations of motion which simultaneously describe the molecular and quantum transport dynamics, and obey the 
fluctuation/dissipation relationship. We also introduce an algorithm which yields the most probable molecular 
and quantum transport pathways in rare, thermally-activated reactions. In the second part of the paper, we apply 
this formalism to simulate the propagation of a quantum charge during the collapse of a polymer from an initial 
stretched conformation to a final globular state. We find that the charge dynamics is quenched when the chain 
reaches a molten globule state. Using random matrix theory we show that this transition is due to an increase of 
quantum localization driven by dynamical disorder. 



I. INTRODUCTION 

Understanding the mechanisms involved in the transport 
of charged and neutral quantum excitations inside macro- 
molecules is a key step towards realizing nano-scale organic 
devices with functional (opto-)electronical properties, notably 
molecular wires and antennas. This perspective has motivated 
a huge activity devoted to investigating the conductance of 
inorganic^] organic^ and biological 8 11 polymers. 

Unlike quantum wires made of solid-state nano rib- 
flexible molecular wires in solution can undergo 
conformational transitions. This feature raises the problem 
of understanding the implications of the conformational dy- 
namics on the quantum transport properties of the molecule. 
Recent experimental and theoretical studies have shown that 
even small thermal fluctuations of the DNA backbone can sig- 
nificantly alter its conductivity ED These effects should be 
greatly amplified in more flexible polymers, which can un- 
dergo cooperative transitions and major re-arrangements of 
the three-dimensional structure. 

Investigating the quantum charge transport dynamics dur- 
ing a coil-globule transition of a flexible chain provides an 
ideal framework to probe the conformation/conductance rela- 
tionship and assess the role of quenched and dynamical disor- 
der. In particular, studying such a process involves the explo- 
ration of the crossover from a regime in which the quantum 
transport dynamics is effectively one-dimensional to one in 
which quantum excitations can diffuse in three dimensions. In 
addition, understanding under which conditions a molecular 
wire can display different conducting behavior in the swollen 
and collapsed phases may lead to designing nano-scale molec- 
ular switches, which can be activated by means of chemically- 
or thermally- induced unfolding. 

From a theoretical and microscopic standpoint, the study of 
quantum transport in molecular systems requires a formalism 
in which the electronic excitations are explicitly taken into 
account and the atomi c nu clei are coupled to a solvent. In 
the existing approaches^E^EDD the equations of motion which 
describe the time-evolution of the charge density and atomic 
coordinates are postulated phenomenologically. For example, 
in Ref. 11 the DNA conductivity was computed by assuming 
that the atomic nuclei obey Newton's equation with a clas- 



sical inter-atomic force field. In addition, charged groups 
in the DNA backbone and counter ions in the solvent were 
coupled electrostatically with the quantum charge. Such an 
approach neglects all non-Coulombic interactions between 
quantum and classical degrees of freedom. For example, the 
molecule may lower the total energy by assuming configura- 
tions in which the ionization energy of the quantum charges 
is increased. This interaction has in principle implications on 
the molecular dynamics. 

In view of these considerations, it would be valuable if 
the equations describing the dynamics of electronic excita- 
tions and nuclear coordinates were rigorously derived starting 
from the general theoretical framework of open quantum sys- 
tems^. On the one hand, this would guarantee that all the cor- 
relations between the quantum and classical degrees of free- 
dom are consistently taken into account. On the other hand, 
it would ensure the fluctuation-dissipation relationship is re- 
spected, hence that the correct thermodynamics is recovered 
in the long-time limit. 

In the first part of this paper, we use the Feynman- Vernon 
path integral formalism to provide such a derivation and de- 
velop a rigorous microscopic theory for the dynamics of 
quantum excitations in molecular systems in solution. The 
same path integral formalism is also used to derive an algo- 
rithm which efficiently yields the most probable molecular 
and quantum transport pathways in rare thermally-activated 
reactions. 

In the second part of this paper, we use our quantum and 
stochastic equations of motion to investigate the propagation 
of a quantum charge inside a homo-polymer undergoing a 
coil-globule transition. We find that the charge dynamics is 
strongly suppressed due to a quantum localization driven by 
the dynamical disorder. Interestingly, this effect sets in only 
when the polymer reaches a compact conformation. 

The paper is organized as follows. In section [II] we de- 
fine define the quantum Hamiltonian for system. In section, 
[III] we construct the path integral representation of the time- 
dependent probability to simultaneously observe the molecule 
at a given conformation and the charge at a given molecu- 
lar site. This path integral is then analyzed in section |TV] in 
saddle-point approximation, and in section [V] we derive the 
set of equations of motion for the system and an algorithm 
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to predict the most probable pathways in thermally activated 
transitions. Section[Vl]is devoted to the study of the collapse 
of the chain. Results, conclusions and outlooks are summa- 
rized in section [VHl 



II. A MICROSCOPIC MODEL FOR QUANTUM 
TRANSPORT IN MACROMOLECULES 

Quantum transport processes in dynamical molecular sys- 
tems depend on the time-dependent structure of the electronic 
ground-state and excited states. Clearly, an approach in which 
these wave-functions and the nuclear coordinates are self- 
consistently calculated ab-initio would be extremely compu- 
tationally expensive. 

A commonly efficient strategy to reduce this computational 
complexity consists in coarse-graining of the electronic prob- 
lem into an effective configuration-dependent tight-binding 
Hamiltonian. In particular, in the so-called Fragment Or- 
bital approach^, one first identifies a set of molecular sites 
(so-called fragments) where the charge (i.e. a hole) can be 
found. For example, in DNA these sites can be identified with 
with individual Watson-Crick bases or with base-pairs. Then, 
the highest occupied molecular orbitals for each fragment — 
herby denoted with (j), (/ = 1, . . .N s ) — are computed in den- 
sity functional theory and Born-Oppenheimer approximation, 
neglecting the effect of the rest of the molecule. 

The hopping of the charge between different molecular 
fragments is controlled by the standard tight-binding Hamil- 
tonian: 



Hmc = E fim &i a 



(1) 
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tained from the fragment orbitals |(j)/) and |0,„): 
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where Hks is the Kohn-Shawn Hamiltonian. The hopping pa- 
rameters Tj m and the on-site energies e/ can be taken to be 
real-valued. 

In addition, if the molecule is connected to electrodes, the 
Hamiltonian Hmc must include also the coupling with the 
donor and the acceptor: 



Hmc -> H MC + {T m a\a D + h.c.) + {T NsA a\a Ns +h.c) 



(4) 



In the Born-Oppenheimer approximation, the fragment or- 
bitals in Eq.s Q and ([3]) depend parametrically on the coor- 
dinates of the atomic nuclei, which evolve in time under the 
effect of the inter-atomic forces, and of the interactions with 
quantum charge and with the solvent. 

A consistent way to describe this dynamics is to consider 
the fully quantum Hamiltonian: 



In this equation, Hmc is the tight-binding Hamiltonian defined 
in Eq. while the Hamiltonian Hm controls the conforma- 
tional dynamics of the molecule in the absence of quantum 
excitations, 



a=l 



2M 



-HQ), 



(6) 



where Q — (qi,. . . , q™ ) is the set of all N p atomic coordinates 
and V (Q) is the molecular potential energy which includes the 
interaction between the different atoms and possibly a term to 
account for the electrostatic and hydro-phobic/philic interac- 
tion with the solvent. 

The part of the Hamiltonian Hg + Hmb describes the cou- 
pling of the molecule with a thermal heat-bath, modeled with 
an infinite set of harmonic-oscillators: 

N " - /ft 2 1 \ 

N p oo / c 2 \ 

Hmb = E E -cfiifa + ^-7a& ■ ^ 

X = (xi,X2, ■ ■ ■) and II = (711,712, • • ■) are the harmonic oscil- 
lator coordinates and momenta, /jj and co ; - denote their masses 
and frequencies and ci are the couplings between atomic and 
heat -bath variables. The last term in Eq. ([8j is a standard 
counter-term introduced to compensate the renormalization of 
the molecular potential energy which occurs when the heat- 
bath variables are traced out (see e.g. discussion in Ref! 2 ^!). 

The Hamiltonian |5]) describes a close system at the fully 
quantum level. In the next sections, we shall use the path 
integral formalism to trace out the heat -bath variables and take 
the classical limit for the nuclear degrees of freedom. 



III. PATH INTEGRAL REPRESENTATION OF THE 
QUANTUM-DIFFUSIVE DYNAMICS OF THE SYSTEM 

In this section we derive a the path integral representation 
of the time-evolution of the system described by the Hamilto- 
nian (j5J. 

Let us assume that the molecule is prepared in some config- 
uration Qo and that a hole is initially created at some monomer 
site Icq. We are interested in computing the conditional proba- 
bility Pt{kf,Qf\ko,Qo) mat after a time interval t the molecule 
is found in conformation Qf and the charge at the site kf. Such 
a probability is described by the following time-dependent re- 
duced quantum density matrix: 



P t (kf,Q f ,\k 0l Q ) - 



Tr[\k f: Q f )(k f Q f \p(t)} 
Trp(f) 

Tr[\k f ,Q f )(Q f ,k f \e- 



p(0) eh 



Hn 



Trp(O) 



(9) 



where p(0) = \Qoko) {Qoko\ e k b tI1b is the initial density ma- 
trix, which assumes factorization with a thermal distribution 
for the heat-bath variables. 



H = Hmc + Hm + Hb+ Hmb ■ 



(5) 
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The goal of this section is to represent the conditional prob- 
ability |9]l as a path integral. This can be done relying on the 
so-called Feynman- Vernon formalism for open quantum sys- 
tems, which has been extensively applied to study quantum 
brownian motion 2 ^. 

The path integral representation of the density matrix can be 
constructed in a way which is conceptually similar to that used 
to represent the standard Feynman propagator. The main dif- 
ference is that it requires to perform the Trotter decomposition 
of two real-time evolution operators (forward and backwards 
in time) and an one imaginary time evolution (initial thermal 
distribution of the heat-bath variables). 

It is convenient to adopt a field-theoretic language to repre- 
sent the dynamics of the quantum charge, while using the po- 
sition representation to describe the evolution of the atoms in 
the molecule and of the harmonic oscillators in the heat-bath. 
In practice, this corresponds to introducing the following reso- 
lution of the identity at each slide in the Trotter decomposition 
of the real- and imaginary- time evolution operators entering 
Eq. ([9): 



pearing at the exponents read 

~MQ{t') 2 
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The path integrals over the harmonic oscillator variables are 
Gaussian and can be carried out analytically. One obtains 



l dQ l 



f dx I 
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(10) 



where <I> = ((j)i , . . . §n s ) are the eigenvalues of the bosonic co- 
herent states constructed from creation and annihilation oper- 
ators in Eq. 

An advantage of adopting a field-theoretic representation 
for the quantum charge dynamics is that the statistical weight 
of the coherent field configurations in the path integral can be 

written as the exponent of an action, i.e. / (0(^(0^* eTi^* ,<s> \ 
This property is very useful to develop the saddle-point ap- 
proximation, and is not satisfied if the dynamics of a tight- 
binding Hamiltonian is represented using the discrete position 
eigenstates. 

The conditional probability (|9]l in path integral form reads: 



P t (kf,Q f \ko,Q ) = fJvQ J ^(0^(0) 

e -E m 0,';(o)<t>L(o)+^[e'VVi [ Qf d@' [ 

jQn J 



-ty" k *(Q)£(t) ■ e -L„C(O0;;,(')-K s Mc[e"/.0"i e -4>FviQ',Q"}_ 

(15) 

&fv[Q',Q] is the so-called Feynman- Vernon influence 
functional 25 , which describes the fluctuation and dissipation 
induced by the coupling with the heat-bath and reads: 

<Pfv[Q',Q"} = {j/t' £dt" [Q'(t')-Q"(t')) 
[L(t'-t")Q'(t")-L*(t'-t")Q"(t")] 



Q' 2 (t')-Q"\t') 



m.col 



(16) 



P t (k f ,Q f \ko,Qo) = JdX J dX { J dX 2 'DXe'^ 
<DX' f\x" f Qf <DQ' [ V4v4* *L(0^(0) 



e -E m *L*(0)^(0) e i{SMclQ ■* .<f *]+Smb[Q'.x'}) 
J <DQ"J <D$'<D$ >,;(0^ (0) 



(11) 



In this equation <|)j (t) are the complex bosonic fields as- 
sociated to the charge coherent states propagating forward and 
backwards in time respectively and the action functionals ap- 



L{t) is a two-point correlation function which encodes the 
physics of the coupling of the molecular coordinates with the 
heat-bath and reads: 



C 2 : 



coth 



I cosfco.f ) — ; sinfco.f 1 

2k R T I y J ' y 1 



.(17) 



Note that the strength of the fluctuation and dissipation in- 
duced by the solvent and the time scales at which memory 
effects die out can be tuned by changing the parameters in the 
harmonic bath Hamiltonian (|8j. In particular, here we con- 
sider the so-called ohmic bath limit ( see e.g. RefP3) , in which 
the L(t) reduces to 



L{t) -> L ohm (t) 



-n^ m+ ^jt m (18) 
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and y is interpreted as the friction coefficient. In section|V]we where the functionals Q kf and M kf read 
shall show that this choice leads to the natural generalization 
of the classical over-damped Langevin dynamics. 



IV. SADDLE-POINT APPROXIMATION OF THE PATH 
INTEGRAL 

The path integral ( [T5] > provides an exact representation of 
the conditional probability given in Eq. Q and cannot be 
solved without relying on some approximation. To this end, 
we observe that for any time t there exists an obvious sum 
rule: 



Qk f [o] 
M kf [a"} 



-E m «(0)<t>,„(0) 



(25) 



" II ^ If j. 



c /,-Vv,/ .0 o ] 

while the functional 'W and Sqm are defined as 



(26) 



W[x,y,p ,p ] = 



J^t'{Mry-v(r + l) + v( r -l) 



£ / dQ f P t (k f ,Q f \ko,Qo) = L 



(19) 



i, m 



fin 



9lm-flm[r-- \ P, m 



The idea is then to use the path integral representation ( 15 i to 
implement this condition, and analyze it in saddle-point ap- 
proximation. 

We begin by changing variables for the molecular coordi- 
nates: 



E / *W) 



ih ~g§i m -Gim{t') 



(27) 

MO- 

(28) 



The path integral ( 24 1 is still and exact representation of the 
(20) sum-rule. The saddle-point approximation is implemented by 
imposing the stationarity of the exponents with respect to the 
Next, we use the path integral ( fTB) to re-write the sum-rule tensor fields a lm , a lm , p lm and p /m and with respect to the 



y (0 = e'(f)-e"(0 r{t) = \{d{t)+d'{t)) 



1 = E / Drfdyf 1>4d4* ^(f)^(0) 
k f JQn JO J 

(21) 

where the new expression for the Feynman- Vernon functional 
reads 



molecular paths y and r. In particular: 

Imposing the stationarity with respect to the r path leads 
to the equation 



My = 2Myy-2, 



V r 



y 



-Vir- 



y 



l.m 



dr L 



y 



flm(r+~) Pi m -fim(r- 



y 



P/m 



(29) 



<P FV [r,y}-j Q d t '[-^y^) + ^r.y 



iMj 



• Imposing the stationarity with respect to the density ten- 
sor fields p lm and p jm leads to the equations 



(22) 



We now introduce a set of tensor fields P; OT (0>Pzm(0> a / m ( f ) 
and a', (t) into the path integral by means the functional iden- 



Glm — flm 
°lm = flm 



(30) 
(31) 



titles 



1 = j (Da (Dp e * E,J " S '° d '' c ' lm ( p '<» ' ) , 
1 = f W'©p% _ ^ E '' m/ ° rfr ' "'"( p ''" _ *''' <t> '" ). 



Imposing the stationarity with respect to the conjugate 
fields a, (t 1 ) and a, (t 1 ) leads to the equations 



pL(0 = 



1 



LkfM^a"} Q k i[a'] 



E M ^/[ 



(23) 



The new expression for the sum-rule ( [T9| l is then 

" <Dr J rDye^wAe-^rvM [^%\a] M kf [a'] ](24) 



(32) 



pL(0 



77 E / 



♦L*wt?"(ot:(oti(o)«-*** f » v *' a * 1 



iQo 



(33) 
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The set of saddle-point equations of motion (29 i-(|32]i are 
simultaneously satisfied if, for any t' £ [Q,t] one imposes 



y(0 = o, 

<4i(0 = °'L( t ')=flm[r(t% 

p'to,(0 = pL(0 = pi«(0- 



(34) 
(35) 
(36) 



Hence, backwards and forward evolution coincide, at the 
saddle-point level. 



Using Eq. ( 32 1 we note that p/„, can be re-written as: 



Plm(t') = (¥r a i| a\a m \Vu 



(37) 
(38) 



where the quantum state is defined as: 

|¥ [a] ) = (39) 
and the time-dependent Hamiltonian H e ff[a] is defined as 

I .111 

Hence, the field p/ m is identified with the (reduced) density 
matrix, evaluated on the state obtained by evolving for a 
time t' the initial quantum state |&q), according to the time- 
dependent Hamiltonian H e ff[a}. 



QUANTUM AND STOCHASTIC EQUATIONS OF 
MOTION 



In this section, we use the path integral representation of the 
conditional probability (|9]l and the saddle-point relationships 
(|34|>-([36ll to derive a set of equations which describe the evo- 
lution of the quantum charge and the classical atomic nuclei, 
for a molecule in solution. 

Our strategy consists in estimating the path integral over the 
charge density field p and its conjugate a field in the lowest- 
order saddle-point approximation developed in the previous 
section. On the other hand, the integral over the molecular 
coordinates y is evaluated at the one-loop level, i.e. includ- 
ing the effects of leading-order fluctuations around the saddle- 
point solution y(t) — 0. This guarantees that the dynamics of 
the molecular coordinates r is stochastic, even at the classical 
level. 

To implement this program, we impose that the density ma- 



trix field p /m = p /m = p/,„ and its conjugate field a 



lm 



0/,„ obey the saddle-point relationship ( |35) , po} and we focus 
on the remaining part of the path integral, which concerns the 
molecular degrees of freedom, T>r and 

*(Qf\QoM= Porftye^W**™), (41) 
JQo J 



where we have used Eq. ( 35 1 to eliminate a; m . 



Retaining only the lowest orders in the expansion of the 
y(f') path around the saddle-point solution y(t) = the func- 
tional at the exponent reads 

W[r,y,p]+ih $> FV [r,y\ = 
= -J dt' [y{t') ■ (Mr(t')+Myr(t') 



+ -< I /[r(0,p(0]) + 



i k B TMy 



y 2 (t) 



(42) 



where 



V[r,p]=V{r)+Y,Pimflm(r)=V(r)+Tr\f(r)P] (43) 

lm 

Since we have kept only terms which are at most quadratic 
in y(t'), the path integral over this variable can be evaluated 
analytically by completing the square. The resulting (unnor- 
malized) expression for the path integral (41 1 is 



ajG2/|Go;[p]) = 



Q f Sp dt {Mf(t')+Myt(t')+l r nr,P]} 2 

Qo 



(44) 



We emphasize that the action in the exponent generalizes the 
well-known Onsager-Machlup functional 31 which appears in 
the path integral representation of the classical Langevin dy- 
namics (a brief review is given in the appendix). 

In macro-molecular systems in solution inertial effects are 
damped at a time scale 10~ 13 s, which much smaller than the 
time scale associated to local conformational changes. If the 
acceleration term Mr at the exponent is neglected the path in- 



tegral (44 1 becomes 



Qf -A f'dt 



Tire 



My di 



00 



(45) 



which is completely equivalent to the stochastic path integral 
defined in Eq. ( |A1 1[ ) for classical over-damped Langevin dy- 
namics . 

Hence, we conclude that the evolution of the system can 
be described by the following set of quantum and stochastic 
differential equations: 



d r« = (V(r) +Tr[p /»]) +r, a (0 



dt 



(46) 



i[/M,P], 



and we have condensed the Eq.s ( f3"7| ), ( |40] >, ( |35| ) and ( |36| ) into 
a single Van-Neumann Equation, where [p(f)]/m = P/m(0 an d 
[f]im( r ) — fim( r )- T 1 CX (0 is tne usual white delta-correlated 
Gaussian noise of the Langevin dynamics, 



(r,«(r).r,P(0)> =6^ 8«P 5(f). 



(47) 



The set of equations ( 46 1 represents main result of this pa- 
per, as far as the formalism is concerned. Some comment on 
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these equations are in order. First of all, we emphasize that 
the Langevin equation for the molecular coordinates contains 
the force term 



F a (Q) = Tr 



(48) 



This term is of fully quantum origin, as it follows directly 
from the the Hamiltonian Q. It expresses the influence of 
the charge distribution on the molecular motion: the atoms 
are driven towards configurations for which the energy of the 
charge is lower. This type of non-Coulombic charge-nuclei 
interaction is not included in the standard phenomenological 
approaches which have been used to simulate the dynamics of 
molecular wire s I P 6 ! 

If the molecule contains charged atomic groups (e.g. like in 
DNA), one needs also to include their Coulombic interaction 
with the propagating quantum charge. This can be rigorously 
done by extending the potential energy function V(Q) in the 
original Hamiltonain ([6]) to: 

P(G)-> v(G;{«/4}'=wvJ 

i N P N s 
<x=l/=l 

where e is the charge of the hole and q a is the (partial) charge 



<lo.e ,, Q , 
2 ^ ^\ r —r-\ a i a '' < 49 > 

/ a=l/=ll r ' r «l 



of the a-th atom. At the level of the Langevin equation (46 1, 
this coupling produces an additional Coulomb force term in 
the form: 



r c,, s f j e « e P" r 
F c (r,Plm) = \ dr a.~( V2 M '«' 



(50) 



where K/ a = (r a - n)/\r a - r{\. 

The equations of motion ( |46| can be integrated by simulta- 
neously evolving the molecular positions and the initial quan- 
tum state, and then computing the corresponding density ma- 
trix field: 

f r„(f + At) = r„(0 - ^ (V [r(0] +Tr[p/»]) 



(51) 



|*P(? + Af)) =e-% %/[ r Wl|^(?)}, 
p /m (f +Af) = (W(t+At)\l){m\W(t+At)}, 

where Z, a (t) is a stochastic variable sampled from a Gaussian 
distribution with zero average and unitary variance. 

The solution of the set of equations ( |4"6] > through the al- 
gorithm pi) yields the real-time evolution of the atomic nu- 
clei and of the charge density, starting from an initial con- 
dition. Like any approach based on the time integration of 
the equations of motion, this method is expected to be com- 
putationally inefficient to investigate rare thermally-activated 
reactions. The main difficulty arises from the fact that, on 
average, the first reactive event occurs on a time scale which 
scales exponentially with the height of barrier. Hence, very 
long trajectories have to be generated in order to investigate 
the dynamics of the reaction. 



In order to overcome these difficulties, the Dominant Re- 
action Pathways (DRP) method was recently developed 20 24 . 
This approach uses the path integral representation of the 
classical Langevin dynamics to derive a variational princi- 
ple which yields the most probable reaction pathways con- 
necting given initial and final configurations. To date, the 
DRP algorithm has been successfully applied to investigate a 
number of reactions which are extremely hard to simulate us- 
ing standard MD algorithms, including conformational transi- 
tions of peptides from ab-initio simulations^, protein folding 
within realistic atomistic classical models^! as well as chemi- 
cal reactions 29 or particle-surface interactions 30 . Recently, the 
formalism has been extended to include quantum corrections 
to the motion of light nuclei 3 - -! 

The Feynman- Vernon functional integral formulation of the 
dynamics developed in the previous sections offers the proper 
framework to extend the DRP method to the case in which the 
molecule contains quantum excitations. The starting observa- 



tion is that the solutions of the saddle-point equations ( 34 1 



( 36 1 and the molecular path minimizing the effective action in 
the exponent of Eq. ( |44| yield the most probable coupled evo- 
lution of the charge density and of the molecular coordinates 
during transitions. 

Assuming that the initial and final molecular configurations 
and the initial position of the charge are held fixed, the saddle- 
point equations can be solved self-consistently by means of 
the following iterative algorithm: 

1. A reaction pathway r(t) connecting the initial config- 
uration Qo and the final configuration Qf is computed 
by minimizing the Onsager-Machlup action at the ex- 
ponent of Eq. ( p5] l, neglecting the coupling with the 
density matrix pi m (t'). This can be ef ficien tly done by 
using the algorithms described in Ref. d 27 ! 28 -! 

2. The deterministic evolution of the density field p/ m (f') 
is computed from Eq.s ( [39| and ( |40| , where the tight- 
binding coefficients in U e ff are evaluated along the path 
r(t ) obtained in the previous step. 

3. The density matrix obtained in the previous step is 
plugged in Eq. ( |4"5] l, and an improved estimate of the 
molecular reaction path r(t) is obtained by numeri- 
cal relaxation of the generalized Onsager-Machlup ac- 
tion d45ll. 



The process is iterated until convergence. 



VI. CHARGE LOCALIZATION IN A COLLAPSING 
CHAIN 

Having developed the appropriate formalism, we are now 
ready to analyze the quantum transport dynamics during the 
collapse of a homo-polymeric chain. 

According to the standard scaling theory 34 , the existence of 
an Anderson metal-insulator transition (MIT) in the thermo- 
dynamic limit and at zero temperature critically depends on 
the number of spatial dimensions. In particular, in the one- 
dimensional Anderson model, an arbitrarily small amount of 
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disorder is sufficient to stop electric conduction. In three- 
dimensions, a mobility edge is developed and the system re- 
mains in the metallic regime for sufficiently small disorder. If 
such idealized conditions were good approximations for re- 
alistic molecular wires in solution, then one might argue that 
the collapse into a molten globule of a polymer may lead to 
an increase of the conductance. 

On the other hand, physical molecular wires in solution are 
complex mesoscopic systems subject to dynamical disorder 
and finite temperature effects. As a consequence, their con- 
ductance depends on a number of specific physical conditions, 
such as e.g. the chain length, heat-bath temperature, noise 
memory function, and so on. 

Nevertheless, this systems may still display different con- 
ducting properties, depending if the polymer is in a coil con- 
figuration (where the transport dynamics is effectively one- 
dimensional) or in a compact phase ( where the charge can 
diffuse in three-dimensions). 

This discussion raises two main questions: does the dynam- 
ical disorder driven by the conformational fluctuations of the 
chain induce a strong charge localization already in the coil 
state? Is such a localization enhanced or reduced when the 
system reaches a compact configuration? 

To address these questions, we use the framework devel- 
oped in the previous sections to simulate the charge migration 
in a simple coarse-grained model for a polymeric collapse. 



A. Definition of the Model 

We consider a chain composed by N p = 36 beads, each 
one representative of an individual monomer. These particles 
are linked by harmonic springs and can interact at distance 
through Van-der-Waals potential, which simulates the effec- 
tive hydrophobic attraction and the steric repulsion: 

1 



ne) = ^l>(k+i-'al- fl ) 2 

* a 



+ 2l 4£ 
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(52) 



The hole can propagate by hopping across N s sites, which 
are identified with the chain's monomers, hence N s — N p . The 
chain is not coupled to electrodes. 

In the present illustrative toy model, the couplings in the 
effective configuration-dependent tight-binding Hamiltonian 
([l} are not evaluated microscopically from electronic struc- 
ture calculations. Instead, their dependence on the monomer 
coordinatesis is defined phenomenologically as follows: 



Tij(Q) = T e 



ei{Q) = e . 



(53) 



Note that the charge can hop between monomers which are 
specially close, even if they are far in the chain sequence. 

The numerical value of the constants entering Eq.s p2| ) and 
(53 i are listed in table [I] We consider a heat bath of tempera- 
ture T = 300 K and, while the mass and the friction coefficient 



TABLE I: Parameter of the Hamiltonian for the toy-model molecular 
wire 





a 


e 


a 


[kJ mol~' nm~ 2 ] 


[nm] 


[kJmor 1 ] 


[nm] 


1000 


0.38 


4 


0.4 


To 


£o 


aj 




[kJmol" 1 ] 


[kJmol -1 ] 


[nm] 




2 


-200 


0.3 





are chosen in such a way that My = 2000 amu ps~'. This set 
of parameters was chosen to ensure that the end-to-end propa- 
gation of the quantum charge across the wire and the collapse 
of the molecule occur at comparable time scales. 



B. Quantum Transport and Conformational Dynamics 

The potential energy function of this model is characterized 
by a low degree of frustration. As a result, the collapse of the 
chain from an initial stretched coil configuration to a compact 
globular state occurs in just few ns. Such a time interval can be 
simulated directly by integrating the equation of motion ( ftp} . 

Let us consider a polymer which is initially prepared in the 
fully stretched configuration shown in Fig. [1] with a quan- 
tum charge localized at the first monomer = 0)) = |1). 
Then, the evolution of the system is predicted by means of the 
algorithm ( |5T| ). 

The elementary time-step in the integration of the Langevin 
equation for this simple coarse-grained model can be cho- 
sen to be AtMD = 10 fs, while the time-step in the quantum 
evolution of the charge wave-function has to be chosen much 
smaller, e.g. Afg = 1 / lOAf^D- With such a choice of parame- 
ters, simulating few ns of dynamics takes only a few minutes 
on a standard laptop computer and the normalization of the 
charge wave-function remains constant within 1% accuracy 
over the entire simulated time interval. 

Let us begin by discussing the charge migration dynam- 
ics in a static wire, i.e. assuming that the molecule remains 
'frozen' in the initial fully stretched configuration. This con- 
figuration corresponds to a complete absence of disorder. In 
Fig. [T|we report the calculated evolution of the charge density 
in few ns of dynamics. We note that the propagation of the 
charge occurs through an attenuating traveling density wave, 
which reaches the end of the wire in about 3 ns. The sec- 
ondary peaks are due to the fact that, at any instant, the charge 
has finite probability not to perform the transition. In the long 
time limit, the primary peak looses intensity, which signals 
that the charge is delocalized, with a consequent increase of 
the entropy. 

Fig. [2] shows that the quantum transport dynamics is rad- 
ically changed if the polymer is allowed to move under the 
effect of the internal forces and of the fluctuations induced by 
the solvent, which introduce dynamical disorder in the tight- 
binding Hamiltonian. We can clearly see that the migration of 
the charge is strongly hindered and that the initial density peak 
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Monomer index 



Monomer index 



1.8 ns 



2.6 ns 



Monomer index 



Monomer index 



FIG. 1 : Time evolution of the charge density distribution in a frozen 
fully stretched polymer. The pictures in the upper-right corner repre- 
sent the polymer conformation, the color of the backbone is propor- 
tional to the charge density. 













0.2 ns 


* 1 [ I 












Monomer index 





0.8 ns 


J m 




V W ..-v.... 









Monomer index 



1.8 ns 




Monomer index 



Monomer index 



FIG. 2: Time evolution of the charge density distribution in collaps- 
ing polymer. The pictures in the upper-right corner represent the 
polymer conformation, the color of the backbone is proportional to 
the charge density. 



obey the so-called Wigner-Dyson distribution 1 , 



7ZS 

Pwd(s) = —t 



(54) 



In particular, we note that Pwd{s) vanishes for s — > 0. This 
feature is called level repulsion and is a characteristic signa- 
ture of delocalized eigenstates (i.e. of the metallic regime). 

By contrast, the localization of the wave-functions reduces 
the overlap between adjacent modes, hence suppressing the 
level repulsion. In particular, in the extreme case in which 
there is no overlap between neighboring eigenvectors, the 
level-spacing distribution follows a Poisson distribution: 



Pp(s) = e-\ 



(55) 



We have computed the level spacing distribution by simu- 
lating 48 independent 3 ns-long evolutions starting from the 
same stretched coil configuration, with the charge initially lo- 
calized at the left terminus of the chain. The level spacing 
distribution at different instants during the transition was then 
obtained by explicitally diagonalizing the quantum Hamilto- 
nian ([l}. 

The results are summarized in Fig. [3] In the early time 
of the collapse the spectrum is well explained by a Wigner- 
Dyson statistics, hence the system displays good conducting 
properties even in the presence of disorder. On the other hand, 
after about 1 ns, revel repellence has completely disappeared 
and the spectrum is well described by a Poisson distribution, 
hence the system behaves like an insulator. 

Note that the observed behavior contrasts with what would 
be expected from scaling arguments, implying that the meso- 
scopic character of this system and thermal effects play an 
important role. 

The recent developments in single-molecule and ensem- 
ble manipulation techniques may allow to look for this ef- 
fect. For example, using a conducting atomic force mi- 
croscopy it is in principle possible to simultane ously measure 
the conductivity^^! and the end-point distance^ED f single 
polymers. Alternatively, by simultaneously applying laser- 
induced temperature jumps 42 and electron pulse radiolysis 43 
on an ensemble of neutral collapsed polymers it is possible 
to drive them into non-equilibrium swollen and electrically 
charged states. Then, the evolution of the conductivity in 
these systems during the subsequent thermal relaxation pro- 
cess may be monitored using time-resolved microwave con- 
ductivity techniques^. 



is completely dispersed after about one ns and hardly propa- 
gates into the second half of the chain. Hence, the charge 
transmission properties of the system are strongly affected by 
the configurational dynamics. 

These results raise the question if the disorder-driven 
charge localization takes place already in the coil state, or 
only when the chain reaches the compact state. This qu estion 
can be addressed in the context of random matrix theorjISEH. 
This theory predicts that the level spacing distribution of a real 
orthogonal Hamiltonian with delocalized eigenstates should 



VII. CONCLUSIONS 

In this paper, we have presented a formalism to quanti- 
tatively investigate the dynamics of quantum excitations in- 
side macromolecules which move under the effect of the in- 
tramolecular forces and of the collisions with the solvent 



1 Note that, as usual, this distribution is normalized in such a way to enforce 
a unit average level spacing. 
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FIG. 3: Level spacing distribution for the tight-binding Hamiltonian, 
evaluated after after 0. 1 ns (upper panel) and after 1 ns (lower panel). . 



molecules. The resulting set of coupled stochastic and quan- 
tum equations of motion ( |46| contains a non-Coulombic force 
term which is not present in the approaches where equations 
of motion are postulated phenomenologically. Such a term 
has a quantum origin and explicitly couples the density ma- 
trix of the quantum particle to the molecular coordinates. The 
numerical simulations performed in a simple coarse-grained 
model have shown that this force can give significant contri- 
bution, of the same order of standard (i.e. Van-der-Waals) 
non-bonded interactions. 

We have also used the path integral formalism for open 
quantum systems to derive an algorithm which yields the most 
probable pathways in the evolution of this system, from a 
given initial to a given final molecular conformation. This 
formalism may provide a computationally efficient method to 
investigate the non-equilibrium quantum transport dynamics 
during rare thermally-activated conformational transitions. 

We have applied this formalism to study the migration of a 
quantum charge, during the coil-globule transition of a poly- 
mer. These calculations have shown that the dynamics of the 
quantum charge is strongly quenched during the collapse. By 
analyzing the statics of the eigenvalues spectrum, we have 
concluded that the disorder-driven increase of localization sets 
in only when the system reaches a compact conformation. 

The formalism developed in this work does not only con- 



cern molecular quantum wires, but can be applied to all 
molecular systems in solution which can support the propa- 
gation of quantum excitations. In particular, in the future it 
would be interesting to use it to investigate the transfer of neu- 
tral atomic excitations between amino-acids in peptide chains. 
This would provide a solid theoretical ground to bridge the 
gap between molecular simulations and experimentally ob- 
servable circular dichroism and Forster Resonance Energy 
Transfer spectra in protein folding. 



Appendix A: Path Integral representation of the classical 
dynamics of a molecule in a thermal bath 

The main goal of this work is to derive a path-integral based 
framework to study the evolution of quantum excitations in 
conformationally evolving molecular systems. 

For comparison purposes, it is instructive to review the cor- 
responding path integral formulation of the classical Langevin 
dynamics which is often used to study the conformational dy- 
namics in the absence of quantum excitations. The so-called 
over-damped Langevin equation reads^* 



qa - 



\ a_ 

Mjdq a 



V (Q) +ii(0, (a = !,...,#,) (Al) 



where y is the friction coefficient, V (Q) is the potential energy 
function entering in Eq. (|6| and r|(f ) is delta-correlated Gaus- 
sian noise, satisfying the fluctuation-dissipation relationship: 



(Tf(0- n p W> 



6k B T 
My 



5 ap S(r-0 (a,p = l,...,A^A2) 



Note that in the original Langevin Eq. there is a mass term 
Mq. However, for macro-molecular systems this term is 
damped at a time scale 10~ 13 s, which much smaller than the 
time scale associated to local conformational changes. 

The stochastic differential Eq. ( |A1| > generates a time- 
dependent probability distribution P(Q,t) which obeys the 
well-known Smoluchowski Eq.: 



dt 



P(Q,t) 



My 



l 

kn~T 



W(Q)P(Q,t) 



By performing the formal substitution 

P(Q,t)-- 



(A3) 



(A4) 



the Smoluchowski Eq. ( |A3[ ) can be recast in the form of an 
imaginary time Schrodinger Eq.: 



HeffV(Q,t) 



where 



H eff = 



My 



v 2 +y e// (e), 



(A5) 



(A6) 



is an effective Hamiltonian operator and 
1 



Veff(Q) = 



4k B TMy 



{(W(Q)) 2 -2k B TV 2 V(Q)). (A7) 
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The conditional probability P(Qf,t\Qo) to find the system 
at the configuration Qf at time t, provided it was prepared 
in the configuration Qq at time t = is the Green's function 
of the Smoluchowski Eq., and can be related to the imaginary 



time propagator of the effective "quantum" Hamiltonian ( A6 1: 

PtiQAQo) = e -^ iV{Q) - ViQo)) (Q f \e-^ff\Q ). (AS) 

Using such a connection, it is immediate to obtain an expres- 
sion of the conditional probability ( |A8| > in the form of a Feyn- 
man path integral 



Pt{Qf\Qo)=e 2k s 



v(Q f )-v(Qo) 

IkaT 



Qj 



where 



V/M = f'dt' 
Jo 



My 



4k B T 



■Veff[r] 



(A9) 



(A10) 



is called the effective action. 

The conditional probability Pt(Qf\Qo) is sometimes written 
also in the following equivalent form: 



Pt(Qf\Qo) 



Qr 
Go 



(All) 



where Som[>"] is the so-called the Onsager-Machlup func- 
tional, 



SomV) = f dt 
Jo 



' . My ( 1 

r+— VV(r) , (A12) 



4k B T V My 



Proving the equivalence between the expressions ( A9 1 and 



(All I is not straightforward, since it involves elements of 
stochastic calculus 45 46 . 
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